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Abstract. The order of a phase transition is usually determined by the nature of 
the symmetry breaking at the phase transition point and the dimension of the model 
under consideration. For instance, g-state Potts models in two dimensions display a 
second order, continuous transition for q = 2, 3, 4 and first order for higher q. 

Tamura et al recently introduced Potts models with "invisible" states which 
contribute to the entropy but not the internal energy and noted that adding such 
invisible states could transmute continuous transitions into first order transitions [l}|4j . 
This was observed both in a Bragg- Williams type mean-field calculation and 2D Monte- 
Carlo simulations. It was suggested that the invisible state mechanism for transmuting 
the order of a transition might play a role where transition orders inconsistent with 
the usual scheme had been observed. 

In this paper we note that an alternative mean-field approach employing 3-regular 
random ( "thin" ) graphs also displays this change in the order of the transition as the 
number of invisible states is varied, although the number of states required to effect the 
transmutation, 17 invisible states when there are 2 visible states, is much higher than 
in the Bragg- Williams case. The calculation proceeds by using the equivalence of the 
Potts model with 2 visible and r invisible states to the Blume-Emery-Griffiths (BEG) 
model, so a by-product is the solution of the BEG model on thin random graphs. 



Potts Models with (17) Invisible States on Thin Graphs 



2 



1. Introduction 

The idea that spontaneous breaking of symmetries gives rise to phase transitions plays 
an important role in our understanding of phase transitions in statistical physics, particle 
physics and even computer science. A canonical model for understanding the breaking 
of g-fold symmetries is the g-state Potts model (5), whose Hamiltonian for uniform 
ferromagnetic couplings may be written as 

= (i) 

(ij) 

where the interactions are along the edges (ij) of some graph, the spins Oi at each vertex 
may take q values and 5 is the Kronecker delta. With this Hamiltonian like spins interact 
with one strength and unlike spins with a different strength. As is well known the phase 
transition for the Potts model in two dimensions is continuous, i.e. second order, if 
q = 2, 3, 4 and first order for larger q. It is generally expected that a phase transition 
associated with the spontaneous breaking of a g-fold symmetry in two dimensions will 
have the same order as the ferromagnetic g-state Potts model, but recent studies have 
suggested this is by no means always the case (6j-j8j. 

It was observed in [l|^4] that one possible mechanism giving rise to such non- 
canonical orders for phase transitions might be the presence of "invisible" states which 
did not contribute to the internal energy but still contributed to the entropy. It is 

possible to write down a Potts Hamiltonian that incorporates such states as 

<? 

^(?,r) = -£<W, J2S Sita 5 Sj>a , 8i = l,---,q,q+l,---,q + r, (2) 

(id) a=1 

which has q visible states and r invisible states for each spin Sj. In the sequel we will 
denote this as the (g,r)-state Potts model. 

A useful way to rewrite this Hamiltonian is to introduce new spin variables o~i where 
(jj = Si if Si = 1, • • • , g and <7j = otherwise p]. The resulting Hamiltonian 

<? 

^(?,r) = 'E^'°*'i- a - Tlnr E^,o, <J% = 0, 1, (3) 

(id) a=1 i 

contains a temperature dependent field term which gives the entropy contribution of the 
invisible states and the additional Kronecker deltas in the first interaction term ensure 
the invisible states do not contribute to the internal energy. 

In the case of a Potts model with 2 visible states and r invisible states we may 
map this Hamiltonian in turn to that of a spin one model, the Blume-Emery-Griffiths 
(BEG) 1 9] model with equal couplings 

(id) (id) * 

U =+1,0,-1 (4) 

where \x = Tlnr. The generic BEG model admits non-equal couplings for the t{tj and 
tit 2 , terms, but this is unnecessary for our purposes. The phase diagram for the BEG 
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model has been studied by numerous authors using various means. In a Bragg- 
Williams mean-field approximation was used, giving a phase diagram whose general 
form is shown in Fig. (1) where a magnetized phase exists in the small //, T region. 



(b) 






. Wa) 



T 



Figure 1. The schematic drawing of the (equal coupling) mean- field BEG model 
phase diagram in the /i, T plane. The second order transition line is shown in bold 
and the first order dashed. The arrowed sloping lines (a), (b) represent different values 
of r when n — TTnr to make contact with the (2,r)-state Potts model. 



The parameter /i is usually called the crystal field as a consequence of the model's 
phenomeno logical roots. 

The mechanism by which the invisible states may affect the order of the transition 
is clear from the diagram. For smaller r as T is increased the system will move along 
a line similar to the the lower of the two sloping lines (since [i = Tlnr), which we 
have denoted by (a) in Fig. (1). It will thus cross the second order portion of the phase 
transition line. For sufficiently large r the line is steep enough to cross the first order 
portion of the transition line, as shown in (b). Quantitatively, a second order transition 
occurs if (q,r) = (2, 1), (2,2), (2,3) and the transition becomes first order for larger r. 
The existence of the low temperature symmetry broken phase and first order transition 
for the (g,r)-state model has also been shown rigorously using random cluster methods 
in [To] for q > 1 and sufficiently large r. 

In this paper we employ another mean-field like model to study the (2,r) state 
Potts model by making use of the fact that spin models on regular random graphs 
(random graphs where each vertex has the same valency) will generically display mean- 



field behaviour since such graphs look locally like a Bethe lattice 11 . Unlike the Bethe 



lattice regular random graphs do have large loops 12 , which renders them closed and 
obviates the need to deal with boundary conditions, as one must for the Bethe lattice 
proper. The weightings to use for spin models living on such graphs can be determined 
by using similar methods to those used for "fat" random graphs, where the perturbative 
expansion of matrix integrals over N x N matrices in the limit N — > oo give planar 
random graphs decorated with the appropriate Boltzmann weights. The regular random 



graphs considered here are generated in the N — > 1, scalar limit of such integrals 13 14 
so we denote them as "thin" random graphs. 
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2. Thin graphs 

The required ensemble of random graphs can be generated from the Feynman diagram 
expansion of a scalar integral. For instance, if we are simply interested in calculating 
the number of undecorated 3-regular random graphs with n vertices we could evaluate 
the integral 



which generates the graphs in a perturbative expansion with a unit propagator from the 
1 cj) 2 term tying together the cubic vertices from the | 3 term. The contour integral in 
A picks out the desired size of graph. Expanding the cubic terms in the exponential and 
evaluating the Gaussian integrals gives the well known result 

" W (2'i)! ' 
for the number of 3-regular random graphs with n vertices. 

When the graphs are decorated with spins we introduce additional variables in 
the integral to give the desired weights to edge interactions. For the Ising model, for 
example, the Hamiltonian on any graph is 

Hi = - ( 7 ) 

<ij> 

where the sum is over the edges connecting nearest neighbour sites. The partition 
function 

M V / 

for the Ising model on an ensemble of thin graphs with n vertices is then given by the 



integral 13 



1 / dX I U( ;_|_u< i- 



where the propagator K is given by 

K-J - ( I 7 ) (10) 

with c = exp(— 2/3) and the integrand, which we shall call an "action" since we are 
thinking in terms of Feynman diagrams, is 

Sj = ij^^Vk - ^(0+ + (11) 

a,b 

Inverting the propagator gives the correct Ising interactions between the two spin states, 
which are proxied by the two variables + , 0_ in the integral. The approach extends 
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naturally to other spin models, such as the q-state Potts models where we may take an 
action of the form [15] 

^4e^- c Em-^e^ ( i2 ) 



with 



i<j i=l 

1 



" exp(2/3) + g-2 (13) 
and calculate the equivalent of equ. (j9]) with the q variables fa. The Potts transition is 
found to be first order for q > 3, as one might expect for a mean-field model. 

To evaluate such spin model partition functions on thin graphs in the 
thermodynamic limit we can employ saddle point methods, in which case the vertex 
coupling A may be scaled out and we obtain the leading term by simply solving the 
saddle point equations for the action: dS/dfa = 0. Phase transitions occur when the 
free energies, given to leading order by S itself, of the different solution branches cross 
over as the temperature or other parameters are varied. Using saddle point calculations 
for the thin graphs gives a straightforward approach to mean-field theory for spin models 
in general and it has been possible to show that the saddle point equations for the 
thin graph actions may be transformed into the fixed point equations of the recursion 



relations 16 used to solve various spin models on Bethe lattices. 



3. BEG on thin graphs 

The appropriate S for the BEG model on thin (or indeed fat) graphs is straightforward 



to write down 17, 18 , even for the more general case of a BEG model with unequal 



couplings where the Hamiltonian is given by 

Hbeg = - f E -fE - » E ( l - *i) • ( 14 ) 

Since this is a spin one model, we employ three variables in Sbeg 
Sbeg = \{4>\ + 4>l + 4>l) - a(l - b)(fah + 2 3 ) 

_ &0102 _I(03 + 3 + A0 3 ) (15) 

where the couplings in the action Sbeg are related to those in the Hamiltonian by 

b + a 2 (l -bf 



exp(— (3 J) 



l-a 2 (l-b) 



2 



eM ~ m = [l-a*(i%V+a*(l-6)'] (16) 
a 3 exp(/3/i) = A 

as can been seen by inverting the quadratic terms in Sbeg to obtain the propagator 
and demanding that the edge weights for the different spin configurations match those 



generated by the Hamiltonian of equ. (14). 
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For the purposes of investigating the (2, r) state Potts model it is sufficient to take 
J = K = 1, which means S'beg simplifies to 

Sbeg = \(<& + 02 + 0s) - «(0i03 + 020s) - \{<t>\ + 02 + A 0a) ■ (17) 
Since 6 = when J = K = 1 the relation between the Hamiltonian couplings and those 
in S'beg also simplifies to 

exp(-/3) = 

1 — a 1 

a 3 exp(/3/i) = A . (18) 
Inverting the first of these we obtain the relation between the coupling a and j3 



?TT (19) 

so the physical range of a is given by < a < l/y/2. With the expression for Sbeg in 
hand, the saddle point equations for the model are given by 

<9S B EG , , ,2 n 

— - — = 01 - a0 3 - 0! = 



<9S, 



dcj) 



2 



a 3 - (f) 2 2 = (20) 



^^ = 03-a(0i + 02)-A^ = O 

and these may then be used to obtain the phase diagram. 

The propagator in Sbeg has a canonical form with unit coefficients for each of the 
quadratic terms in equ. dl5[l7| ). This means that the 03 vertices must pick up an extra 
a 3 coefficient to give the correct Boltzmann weights for edges with (0, 0) spins. This 
manifests itself in the relation between A and and the crystal field coupling [i in the 
original Hamiltonian in equs.( 17p8 ). An alternative would be to scale 03 — > 4>-z/a to 
shift all this dependence to the propagator. We have belaboured this point because the 
relation between fi (in the BEG spin model Hamiltonian) and A (in the random graph 
action) determines the number of invisible states r via fi = Tlnr in equ. (18) 

a 3 r = A (21) 

when we think of the BEG model as a representation of the (2, r)-state Potts model. 

Perhaps the most direct way to see the phase structure of the model is to evaluate 
the magnetization along appropriate slices of the two-dimensional phase diagram shown 
in Fig. (1). The natural choice of variables for Sbeg is given by A, a, which we can 
translate back to fi, (5 as appropriate. The magnetization M of a saddle point solution 
may be calculated by substituting the saddle point values of 0i,2,3 into 

0?-02 

M = ,3 3 . 22 

0? + 3 + A0| 

which counts the number of +1 vertices minus the number of —1 vertices, normalized 
by the total number of +1, —1 and vertices. The saddle point values of 01,2,3 may also 
be substituted into Sbeg itself to obtain the free energy. 
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The explicit expressions of the saddle point values of 01,2,3 are not particularly 
illuminating: there are two families of non-trivial solutions, one with <pi ^ <p 2 
corresponding to a magnetized state and one with 0i = 02 corresponding to an un- 
magnetized state. There is also a trivial solution with 0i = 02 = 03 = 0. The expression 
for the saddle point action S sp in the case of the magnetized solution(s) is, however, quite 
simple 



J sp 



1 1 

6 + 12A 2 



a 1 , 
2A 12A 2 ' 



1 - 4Aa) 



3/2 



(23) 



and makes it clear that we might expect critical behaviour (at least) along the curve 
4A c a c = 1. The saddle point magnetization also confirms this since it is given along 
the various branches by 



M = ± 



A(A -2a± 2a VI -4Aa)(2A - a ± ay/1 - 4Aa) 



(24) 

2A 2 - 6Aa ± 4AaVl - 4Aa + 1 =f y/1 - 4Aa 

To scan the phase diagram it is convenient to evaluate the magnetization M at fixed 

a as A is varied. In Fig. (2) we have taken a = 0.1 (i.e. a low temperature value) and 

plotted all the non-zero branches of the magnetization against A. The branching values 




Figure 2. The magnetization plotted against A at a = 0.1 (i.e. low temperature) 



at the tips of the prongs are at A c = 2.5, confirming that 4A c a c = 1. As A is reduced at 
fixed a (which also corresponds to the crystal field ji decreasing) the magnetization will 
jump at A c since the internal branches of the prong have a higher free energy, signalling 
a first order transition. 

As a is increased the tips of the prong move in towards the horizontal axis and 
by a = 0.37 they have coalesced as shown in Fig. (3) which gives a picture that is 
typical of a second order transition. We can pinpoint the triple point value a t where the 
prongs coalesce and the transition becomes second order by substituting the equation 
for the first order critical line A c = l/4a c into the saddle point magnetizations to get 
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0.7(1 



Figure 3. The non-zero magnetizations plotted against A at a = 0.37 which lies just 
above the triple point value of a t = l/v° ~ 0.3535 . . . 



the magnetization at a c , giving expressions of the form 



Jl - 8a 2 (1 - 2a 2 ) , s 

M ° =± (4ag-l) < 25) 



which show that M c = at the triple point value, at = l/v8 — 0.3535.. with 
A t = l/y/2 ~ 0.7071... 

On the second order transition line the critical value of A decreases as a is increased, 



as can be seen from the full expression for the non-zero branches of M in equ. (24) above 



A(A - 2a + 2a\/l -4Aa)(2A - a + ay/1 - 4Aa) 

M = - t , (26) 

2A 2 - 6Aa + 4Aa v / l - 4Aa + 1 - yj\ - 4Aa 

where the physical solutions for the second order branches are given by choosing the 



upper signs for the various square roots in equ. (24). As a — > 1/2 the critical value of A 
tends to zero and no magnetized solution exists for the remainder of the physical range 
1/2 < a < l/y/2. 

Remembering that increasing a corresponds to increasing T the general picture 
presented from the thin graph saddle point calculation is thus perfectly consistent with 
that observed from the Bragg- Williams mean-field approach shown in Fig. (1). For 
< a < a t = 1/ \/8 first order transitions are observed as A is varied, whereas between 
l/y/8 < a < 1/2 second order transitions are observed with the critical value of A 
decreasing with increasing a (i.e. increasing temperature). For 1/2 < a < l/y/2 on the 
other hand, there is no transition. 

Having obtained the phase diagram of the BEG model on thin <f) 3 graphs we are 
now in a position to re-interpret it as the phase diagram of a (2,r)-state Potts model on 



the same graphs. The key equation here is equ. (21) which states the relation between 
the number of invisible states and the parameters in Sbeg; namely a 3 r = A. The key 
values to insert in this are those of the triple point where the transition changes from 
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being second to first order, a t = 1/ and A t = 1/ \/2, so we find that r at this point 
has the surprisingly large value of 16. 

In summary, on 3-regular random graphs the (2,r)-state Potts model requires 17 or 
more invisible states to transmute the second order transition displayed by the 2-state 
Potts (Ising) model on such graphs into a first order transition. 



4. The Wajnflasz-Pick and the standard Potts model on thin graphs 

Curiously, another mechanism for changing the order of a phase transition discussed by 
Tanaka and Tamura in |4j, the Wajnflasz-Pick model 19 20], bears a close relation to 



the behaviour of the standard Potts model (with no invisible states) on thin graphs 16 



The Wajnflasz-Pick Hamiltonian is Ising-like, but there are g + positive spin states and 
g_ negative spin states 

H W p = - J ^ SjSj - h 2j Si, Si = +1, ■ ■ ■ , +1 , -1, ■ (27) 

<M> * g+ g. 

This Hamiltonian can be transformed to a standard Ising Hamiltonian in a shifted 
temperature dependent external field |20 



n 



T 



J ^ a i a 3 ~ ( h ~ 7T l0 S ~) X) a » <7i = +l,-l. 



WP 



(28) 



We can see in Fig. (4) that changing the temperature at a fixed external field ho will 




Figure 4. The h, T phase diagram for the Wajnflasz-Pick model. The trajectory of 
the system at fixed external field ho as T is increased is plotted for three ratios of g+ 
to ,g_. 



take the system along a line with no transition, a second order transition (if it goes 
through T c ) or a a first order transition depending on the ratio g+/g~ 
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We now consider the standard q-state Potts model where the Hamiltonian is that 
given in equ. ([TJ 

(ij) 

On thin graphs the Potts models display the mean-field behaviour of a continuous 
transition for q = 2 and first order transitions for q > 2. The model on 3-regular 
thin graphs is described by the action of equ. fll2"| , also repeated here for convenience, 

i=l i<j i=l 

and one finds a high temperature saddle-point solution of the form 0j = 1 — (q— l)c, Vz 
bifurcating to a broken symmetry solution 0! = . . . 0, 3 _ 1 7^ g at c = l/(2g — 1) where, 
as before, c = l/(exp(2/3) + g — 2). The magnetization may be defined as 



m = m&f) (29) 

with a corresponding order parameter given by 

(jmax(m) — 1 , , 

M = . 30 

q - 1 

which is zero in the high temperature paramagnetic phase and tends to one in the 
magnetised low temperature phase. 

The graph of m, shown in in Fig. (5) for a 3-state Potts model, is a skewed pitchfork 
but the phase transition takes place neither at O = (q — 1 — 2y/q — l)/(q — l)(q — 5) 
where the low temperature branches meet, nor at P = l/(2q — 1) where the high 
temperature, symmetric branch intersects the low temperature branch, but rather at a 
value Q = (1 — (q — l) _1 / 3 )/(g — 2) which lies between the two where the high and low 



temperature actions S (i.e. free energies) are equal 15 . There is therefore a jump in 



the magnetization of AM = (q — 2)/(q — 1) at this point between the high temperature 
value and that on the magnetized branch. When q = 2 (the Ising case) the pitchfork is 
symmetric and the points 0,P,Q coincide to give a continuous transition 

At first sight, this looks rather different to the picture presented by the Wajnflasz- 
Pick model, but we can recast the Potts model transition on thin graphs as that of an 
Ising model in a temperature dependent external field fl6) by noting that the symmetry 
breaking pattern is q — > (1, q — 1), so setting <fii = <fi and 02 = 03 = • • • = <fiq — 4> m S q 
to give 

~S q = \{q - 1) [1 - c(q - 2)] 2 - ±{q - 1)0 3 + ^0 2 - ^0 3 - c(q - 1)00. 

(31) 



is sufficient to capture the transition. The action of equ. (31) containing only (0, 0) may 
then be rescaled using 

1 
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m 




0.16 0.18 0.20 0.22 0.24 

c 



Figure 5. The magnetisation m for a 3 state Potts model as calculated from the 
saddle point solutions. Only the portion of the graph close to the transition point is 
shown for clarity. 



to give an action equivalent to that of an Ising model in field 

Si = \(<t> 2 + 2 ) - - y0 3 - ^ 3 (33) 

where the new parameters v and k are related to the original c and q via 

1 



(q- 1)V2(1 - (g-2)c) 3 / 2 



c 2 (q - 1) 

' (34) 



l-(g-2)c ■ 

The external field parameter v in equs. ( 33p4 ) depends on the temperature f3 via c, so 
the temperature driven first order transition of the standard Potts model on thin graphs 
appears in this formulation as the first order field driven transition of the Ising model. 
Substituting the value of c at the first order transition point, Q = (1 — (q— l) -1 / 3 )/^— 2), 
into the first of equ. (34) gives v — 1, the zero field point, so the picture of the transition 
in this formulation is similar to that presented by Fig. (4) where the first order transition 
occurs when the system crosses the first order line in the h, T plane. If we take q = 2, 
then v = 1 for all c and the action becomes that of the standard Ising model, which 
displays a second order transition. 



5. Discussion and Conclusions 



We have shown that it is possible to determine the phase diagram of the BEG model 
with equal couplings on 3-regular (0 3 ) random graphs by considering the saddle point 
equations derived from an action which decorates the graphs with the appropriate 
statistical weights. The mean-field phase diagram thus derived shares the properties 
of other mean-field approaches and displays both first and second order transitions. 
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Using the equivalence of the BEG model and a (2,r)-state Potts model with r 
"invisible" states implies that the position of the triple point in the BEG model can then 
be interpreted as showing that a (2,r)-state Potts model requires 17 or more invisible 
states in order to transmute the second order transition of the Ising model on 3 random 
graphs into a first order transition. Although the number of invisible states r needed 
to force a first order transition is larger than the other examples explored by Tamura 
et al, the position of the triple point which determines the critical value of r for such a 
transmutation is not a universal quantity and would be expected to be lattice dependent. 

As further evidence of such non-universality, we could take the in equ. ( T5fl~7 ) 



to be N x N hermitean matrices and consider the N — > oo limit of the integral rather 
than the N — > 1 limit which has effectively been considered here for the thin graphs. 
In this case we would be evaluating the partition function of the BEG model on an 
ensemble of planar random graphs, or if one preferred the partition function of the BEG 



model coupled to 2D quantum gravity. This was done by Fukazawa et al 17 for the 
equal coupling BEG model on 4 planar graphs (i.e. planar, 4-regular random graphs), 
who found that the triple point coupling values on such graphs were given by numerical 
means as approximately a? t = 0.06111 and A t = 0.8329. In this model the relation 
between the crystal field fi and the couplings a and A is given by a 4 exp(/3/i) = A, so 
taking /i = Tlnr as before we have a 4 r = A and a critical r value of approximately 223. 
On such graphs a (2, r)-state Potts model would display a third order transition (the 
order of the Ising transition on such graphs) for up to 223 invisible states. The general 
effect of coupling spin models to planar random graphs is to "soften" the order of the 
phase transition (the Ising transition becomes third order, for example), so many more 
invisible spins appear to be necessary to overcome this compared with, for instance, a 
square 2D lattice. 

We also noted that the first order transition of the standard g-state Potts model on 
thin graphs could be thought of as a different sort of transmutation of the continuous 
Ising transition on such graphs, analogous to the mapping between the Wajnflasz-Pick 
model and the Ising model in an external field also discussed in j3]. We remark in 
closing that the BEG model has been studied on the Bethe lattice using recursion 
relations in |21|. It would be an interesting exercise to confirm that the saddle point 
equations used here were equivalent to the fixed point of these recursion relations, as is 
the case for other thin graph models which have been compared with their Bethe lattice 
counterparts 
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